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Abstract. For an analysis ol microlensing observational data in case of binary 
gravitational lenses as well as for an interpretation of observations of high mag- 
nification events in multiple images of a lensed quasar it is necessary to calculate 
for a given source the microlensing light curve by a fold caustic. This problem 
comes to the numerical calculation of a singular integral. We formulated the 
sufficient condition of a convergence of the integral sum for this singular integral. 
The strictly approach to the problem of a comparison of model results with the 
unequally sampled observational data consists in calculation of the model light 
curve in equidistant points of the canonical dissection of the integration segment 
and a following interpolation of its values at the moments of observations. 



The high angular resolution observations in optical and infrared wave bands 
are of great importance for the solution of many problems of astronomy and 
astrophysics: the search for close binaries, the measurement of angular diameters 
of stars, the investigation of brightness distributions across the stellar disks, 
the study of a fine spatial structure of active galactic nuclei (AGN). A new 
tool to achieving of a high angular resolution has been realized recently - the 
observations of gravitational microlensing. 

In some microlensing events observed by MACHO, OGLE, and PLANET 
collaborations perhaps there are the binary gravitational lenses. For these lenses 
the character of the gravitational potential creates caustics in the source plane. 
The crossing of a caustic leads to sharp variation of a flux received from the 
lensed source. The character of this variation is depend on the source angular 
structure. Thus the observed microlensing light curve contains an information 
about the source brightness distribution (Albrow et al., 1999; Gaudi and Gould, 
1999; Bogdanov and Cherepashchuk, 2000). 

The high magnification events in fluxes received from multiple images of a 
lensed quasar due to gravitational microlensing by stars in a lens galaxy are 
connected also with the complex caustic system which is created by the total 
gravitational potential (Schneider et al.., 1992; Zakharov, 1997). The observa- 
tions of these events allow to obtain an information about the quasar central 
engine - the accretion disk, which surrounds a super massive black hole in the 
galactic nucleus (Gricgcr et al., 1991; Mineshige and Yonehara, 1999; Agol and 
Krolik, 1999; Bogdanov and Cherepashchuk, 2001). The possible angular reso- 
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lution that can be realized in the caustic crossing microlcnsing observations is 
measured by the value of order or less of microsecond of arc. 

From the general considerations we can believe that the case of microlcnsing 
by a fold caustic is a more probable (Gaudi and Gould, 1999). The angular 
sizes of observed stars of the bulge of Galaxy or the Magellanic Clouds as well 
as the sizes of AGN's accretion disks are very small. Therefore we may neglect 
of a caustic curvature and regard the fold caustic as a straightforward line. 
In this case the picture of microlensing depends only on one-dimensional strip 
brightness distribution B(x) in the direction of x axis that is perpendicular to 
caustic. It is known that B(x) is connected with a two-dimensional brightness 
distribution b(x, y) by the integral equation 

oo oo 

B(x)= J J b(£,y)5(£-x)dtdy, (1) 

— oo — oo 

where 5(x) is the Dirac function. For a circularly symmetric source the bright- 
ness depends only on the radial distance to source center r. In this case B{x) 
is connected with b{r) by the Abel's integral equation 

_ ?** 

J \Jr A — x A 

X 

For real astronomical objects the strip brightness distribution is a non-negative, 
continuous, and smooth enough function. 

If £ is the angular distance of point with the brightness B(£) to center 
of source and x is the angular distance of source center to caustic then the 
coefficient of amplification can be written as (Schneider et aL,1992; Zakharov, 
1997): 

A(x-t) = Ao + -^L=H(x-Q, (3) 

where H(x — £) is the Heaviside step function (H = for negative and H = 1 
for non-negative values of its argument), A and K are certain constants for the 
given caustic crossing. The observed microlensing light curve I(x) as a function 
of distance to caustic x can be expressed by the convolution integral equation: 

oo 

I(x)=A(x)*B(x) = J A(x-QB(£)d£. (4) 

— oo 

We can admit further without a lack of community that in the expression 
(3) A n = and K = 1. Far from a caustic the integrand in the expression (4) 
is small enough and we can consider its to be equal to zero for £ < x\ . Then 
as it follows from equations (3) and (4) the estimation of flux value for a given 
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strip brightness distribution when x — x n comes to the numerical calculation of 
the singular integral 

,(*„> . fm*. ( 5) 

J V x n - ? 

Xl 

It is known that the numerical estimation of a singular integral requires of 
certain precautions. The similar problems were examined by Belotzerkovski and 
Lifanov (1985). In particular it can be shown that an attempt to calculate /(xj) 
in arbitrary spaced points Xi when is given in equidistant points & leads 

to a large error of the result. We formulate below the sufficient condition of 
convergence of the integral sum for the singular integral of our type. 

At first we remind the basic definitions and one important theorem. 



Definition 1 



Let the segment [xi, x n ) is divided by points Xj, i = 1, 2, . . . n into n— 1 equal 
sub-segments of length h, and points = 1, 2, . . . n — 1 are the middle points 
of these sub-segments: £j = (xj + x,+i)/2 . Then they say that the sequences 
of points Xi and & create the canonical dissection of this segment. 



Definition 2 



Let a function f(x) is determined in a set D, and for two arbitrary values 
of its argument xi, x 2 € D is valid the inequality 

|/(X!)-/(X 2 )| < fc|X!-X 2 r, 

where fc and a are positive numbers, and < a < 1. Then they say that /(x) 
satisfies the Holder's condition of power a with the coefficient k in the set D. 



Theorem 1 (The mean value theorem) 



Let functions /(x) and g(x) are both integrable in a segment [a, b], and /(x) 
is bounded m < /(x) < M, and sign of <?(x) is invariable in this segment. Then 
the product of these functions f(x)g(x) is also integrable, and 



b b 



J f( x )g( x ) dx = n J g{x) dx, 



where m < fj, < M 



The proof of this theorem can be found in textbooks on mathematical analysis. 
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Take into consideration these definitions and the mean value theorem we can 
formulated the sufficient condition of convergence of integral sum for the rect- 
angle rule that correspond to the singular integral (5) as the following Theorem. 



Theorem 



Let a function B(x) satisfies the Holder's condition of power a in segment 
[xi,x„], and sequences of points x,,i — 1,2,.. .n and = 1,2,. . .n— 1 create 
the canonical dissection of this segment with the step h. Then 



\/Xn - I 



< o{h a+1 ' 2 ) 



Proof 



The sequence of points Xi divides according to Definition 1 the segment 
[xi,x„] into n — 1 sub-segments of length ft. In each of these sub-segments 
the function B(£) is continuous and bounded, and the sign of function g(g) = 
(x„ — C) _1//2 is invariable. Using Theorem 1 and integrate of the function g(£) 
in each sub-segment we can to write 



n-1 



B{^i)h 



i=l i=l & 

where £* is a certain point in sub-segment of number i . Further we have 



A = 



n-1 

E 

»=i 



^(g)v 1 ^ - 5 [^/x7 T^ x7 - yjx n - x t+ i ] - B(£j)h 



Let we multiply and divide simultaneously the expression that closed in 



brackets by factor (y/x n — xi + ^/x n — Xi+{ ). Then we have 



Xi-\-l Xi 



\/~%n ~i \J %n 



< Vft. 



The last inequality is valid so far = ft, and the denominator of this 

fraction is minimal when i = n — 1. 

The common denominator in the expression for A is also minimal when 
i = n — 1, and its minimal value is equal to y/h/2 according to Definition 1. 
Thus we have restriction 
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n-1 



2B($)y/h/2Vh-B(Zi)h 



n-1 



<]TV2X V2B{&)-B{& . 



1=1 



Take into account Definition 2 we have finally A < o{h a+1 ^ 2 ). The Theorem is 
proved. 



The Holder's condition is a more strong restriction on the behavior of a 
function in comparison with a continuity. If the Holder's condition is valid 
then the function is continuous. The opposite statement is generally wrong. 
It is obvious that for implementation of the Holder's condition when a = 1 
sufficiently to demand of bounded first differential of the function. This demand 
is fulfil for strip brightness distributions B{x) of real astronomical objects. 

Remark 2 

For increase of precision of the integral sum calculation in the same sequence 
of points Xi the number of points & can be multiplied by factor k > 1, where 
k is entire. The new sequence can determined as & = xi + h(i — l/2)/k,i = 
l,2,...k(n — 1) . This case is equivalent to the estimation of the singular 
integral in a new sequence of points Xi — Xi + ih/k, i = 1,2,... k(n — 1). The 
new points Xi are partly coincide with points of the previous sequence Xi and 
partly located between them. It is clear that the new sequences Xi and ^ also 
create the canonical dissection of segment [xi,x n ] . 



The observed values of the flux in microlensing observations usually have 
unequal sampling intervals. Therefore the strictly approach to the problem of 
comparison of model results with the observational data consists in the calcu- 
lation of values of the model microlensing light curve in equidistant points of 
the canonical dissection of the integration segment and following interpolation 
of them at the moments of observations. 

Author is grateful to A.F.Zakharov for the attraction of his attention to this 
problem. 
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